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c« ■ Abstract 

We measured atmospheric antiproton spectra in the energy range 0.2 to 3.4 GeV, 
^ , at sea level and at balloon altitude in the atmospheric depth range 4.5 to 26 g/cm 2 . 

The observed energy spectra, including our previous measurements at mountain 
altitude, were compared with estimated spectra calculated on various assumptions 
regarding the energy distribution of antiprotons that interacted with air nuclei. 
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1 Introduction 



Antiprotons are produced in the atmosphere via interaction between cos- 
mic rays and atmospheric nuclei. These atmospheric antiprotons are carry- 
ing important information about the physical processes of their production 
and propagation in the atmosphere. The production and propagation mecha- 
nisms should be similar to those of galactic antiprotons. Therefore the study 
of atmospheric antiprotons is important for understanding the energy spec- 
trum of galactic antiprotons. In addition atmospheric antiprotons are a major 
background for the galactic antiprotons measured at balloon altitude. This 
background has to be estimated by model calculations. However model cal- 
culations of this background [1,2,3,4,5,6] have not yet been verified by direct 
observation. 

We report here the atmospheric antiproton spectra measured with the BESS 
detector in a kinetic energy region of 0.2 - 3.4 GeV at sea level at the atmo- 
spheric depth 994 g/cm 2 in 1997 (BESS-1997), and at balloon altitude over 
the atmospheric depth range 4.5 - 26 g/cm 2 in 2001 (BESS-2001). The an- 
tiprotons observed at balloon altitude are assured not to be primary, but are 
produced inside the atmosphere, because the vertical geomagnetic cut-off en- 
ergy was 3.4 GeV for protons/antiprotons throughout the balloon flight. The 
measured spectra, including our previous data observed at mountain altitude 
at the atmospheric depth 742 g/cm 2 in 1999 [7], will be compared with two 
model calculations. We will discuss how these calculations provide fits to the 
observed data. 



2 The BESS spectrometer 



The detector for the Balloon-borne Experiment with a Superconducting Spec- 
trometer (BESS) was designed [8,9] and developed [10,11,12,13] as a high- 
resolution magnetic-rigidity spectrometer with a large acceptance to perform 
highly sensitive searches for rare cosmic-ray components, as well as precise 
measurements of various cosmic-ray species [14,15,16]. A uniform magnetic 
field of 1 Tesla is produced by a thin superconducting solenoid [17]. The 
magnetic-rigidity (R = Pc/Ze) of an incoming charged particle is measured 
by a tracking system which consists of a jet-type drift chamber and two inner 
drift chambers inside the magnetic field. The deflection (-R -1 ) is calculated for 
each event by applying a circular fit using up to 28 hit points, each with a spa- 
tial resolution of 200 /im. Upper and lower scintillator hodoscopes [13] provide 
time-of-flight and two independent dE/dx measurements. Time resolution of 
each counter is 55 ps, resulting in a 1/(5 resolution of 1.4%, where (5 is defined 
as the particle velocity normalized by the speed of light. The first-level trigger 
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is provided by a coincidence between the top and the bottom scintillators with 
the threshold set at 1/3 of the pulse height for minimum ionizing particles. 
In the BESS-1997 ground experiment, all the first-level triggered events were 
recorded because the trigger rate was 30 Hz, which is low enough to record all 
the events. In the BESS-2001 balloon-flight experiment, the first-level trigger 
rate was too high to record all the events, and thus a second-level trigger was 
issued when a particle's rigidity was calculated to be negative by the onboard 
computer [11,18] to record negatively charged particles preferentially. In ad- 
dition to this biased trigger mode, one out of every four first-level triggered 
events were recorded to provide an unbiased data sample in the balloon flight. 
The BESS instrument also incorporates a threshold-type Cherenkov counter 
[12], whose radiator was a silica-aerogel with a reflective index of 1.03 in BESS- 
1997 and 1.02 in BESS-2001. Antiprotons are distinguished from and fj,~ 
background by imposing that there be no light output from the Cherenkov 
counter. The rejection factor in this analysis was about 50,000 and 7,000 in 
BESS-1997 and BESS-2001, respectively. 



3 Observations 

The BESS-1997 observations at sea level were carried out at the High Energy 
Accelerator Research Organization (KEK), Tsukuba, Japan (36° 12' N, 140° 
6' E), where the geomagnetic cut-off rigidity was 11.2 GV, during the two 
periods of 6th - 11th May, and 7th - 13th December in 1997. The BESS- 
2001 balloon flight was carried out at Ft. Sumner, New Mexico, USA (34° 49' 
N, 104° 22' W) on 24th September 2001. Throughout the flight, the vertical 
geomagnetic cut-off rigidity was about 4.2 GV, which corresponds to a kinetic 
energy of 3.4 GeV for protons/antiprotons. After the balloon reached a normal 
floating altitude of 37 km, where the atmospheric depth is 4.5 g/cm 2 , it began 
to lose altitude and continued descending for 13 hours before the termination 
of the flight at an atmospheric depth of around 30 g/cm 2 . The atmospheric 
depths throughout the observations are shown in Fig. 1. The mean atmospheric 
depth during sea level and balloon altitude observations was 994 g/cm 2 and 
10.7 g/cm 2 , respectively. 



4 Data analysis 

Analysis was performed in the same way as in the previous BESS experi- 
ments [14,19]. We first selected events without interactions inside the BESS 
detector and with good-quality measurements of the rigidity and velocity. 
Particle identification was performed by requiring that dE/dx be consistent 
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with a singly charged particle, no signal was observed from the silica-aerogel 
Cherenkov counter, and particle mass was calculated to be consistent with 
the antiproton mass. In Fig. 2, events between the two curves were identified 
as antiprotons. The number of detected antiproton candidates in the energy- 
range 0.2 - 3.4 GeV were 25 at sea-level in BESS-1997 and 156 at balloon 
altitude in BESS-2001. 

In order to obtain the absolute flux of antiprotons at the top of the BESS 
instrument, we estimated the event selection efficiency, interaction loss prob- 
ability and energy loss inside the instrument, and background contamination. 
The selection efficiency was obtained by using the recorded data, which consist 
mainly of protons since dE/dx and 1/(3 are the same as those of antiprotons. 
The interaction probability and energy loss inside the instrument, as well as 
its geometrical acceptance, were calculated by Monte Carlo (MC) simulation. 
The MC code was tuned and verified by comparing the simulation with an 
accelerator beam test of the BESS detector [20]. The systematic uncertainty 
of the event selection efficiency and antiproton interaction losses were eval- 
uated to be 5% from this beam test. The background contamination due to 
inefficiency of the silica-aerogel Cherenkov counter was estimated by dividing 
N~ by the rejection factor of the Cherenkov counter, where N~ is the number 
of negatively charged particles after applying all selection cuts except for the 
Cherenkov veto to the entire set of events. The rejection factor was estimated 
using the recorded proton events during the balloon flight. In the energy range 
between 1.9 and 3.4 GeV, the background contamination was found to be 20% 
and 5% in BESS-1997 and BESS-2001, respectively. Between 1.0 and 1.9 GeV, 
the contamination was less than 1% for both observations. Below 1.0 GeV, it 
was negligibly small. 

Fig. 3 shows the resultant antiproton flux observed at sea level in BESS-1997 
and at balloon altitude in BESS-2001. These results are summarized in Tables 
1 and 2. The detectable energy range was limited by the threshold of the 
silica-aerogel Cherenkov counter's ability to reject background in BESS-1997. 
As for the BESS-2001 experiment, the upper limit of the energy range was set 
not to exceed the geomagnetic cut-off energy. Results of the model calculations 
discussed below, as well as our previously published data obtained at mountain 
altitude [7], are also shown in Fig. 3. 



5 Model calculations 

Several calculations of the atmospheric antiproton flux have been published [1,2,3,4,5,6]. 
We made a phenomenological calculation of atmospheric antiproton flux un- 
der the observation conditions, following those model calculations. Referring 
to the transport equation given by Stephens [4] , we solved the following equa- 
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tion: 



dJ p {E,x) _ v 9 



+ 



, J p (E',x) Jn(E',x) 

1 — Cm h Q; 

K ' \ in (E>) \ in (E>) 



Jp{E, x) 
HE) 

dE' 



Here, Jp(E, x) is the differential antiproton flux at the atmospheric depth x 
g/cm 2 . The first term on the right-hand side, Qa(Ea, x, E), represents the 
production rate of antiprotons with energy E by the incident particle A with 
energy Ea- Since the flux of incident particles depends on the atmospheric 
depth x, Qa(Ea, x, E) is a function of x. The second term represents the flux 
change due to the ionization energy loss. The average energy loss per g/cm 2 
is indicated by (dE/dx). The third term represents loss of antiprotons due to 
interactions. A(E) is the total inelastic interaction mean free path (mfp). It 
is described as A(E) = [l/\ ann {E) + l/\ in {E)Y\ where \ ann (E) and \ in {E) 
are mfp's for annihilating and non-annihilating processes of antiprotons pass- 
ing through the atmosphere, respectively. The fourth integral term represents 
the tertiary antiproton production rate by non-annihilative inelastic interac- 
tion. $(£", E') signifies the probability that an antiproton with initial energy 
of E' possesses energy E after a collision. The charge exchange probability 
between antiprotons and antineutrons is described as a, which is taken to be 
1/3. Therefore $(E,E')(1 - a) J(E') P / \ in (E') dE' represents the production 
rate of tertiary antiprotons with energy of E by the parent antiproton with 
an energy between E' and E' + dE 1 . In addition to the tertiary production 
from antiprotons, a contribution from antineutrons (n) is also included in this 
integral term. 

The primary cosmic-ray flux used in this calculation was based on the re- 
sults from the BESS-98 [15] and AMS [21] experiments, the results of which 
are in good agreement. The proton, neutron, and helium fluxes at various 
atmospheric depths were obtained by using the transport equation, and a 
contribution from heavier nuclei such as CNO was included, following Papini 
et al. [22]. The antiproton production spectrum by the cosmic rays was taken 
from Stephens's formulation [4]. 

The mfp's, A(E), X m (E) and X ann (E), are the same as those adopted by 
Stephens [4]. These are shown in Fig. 4. The details of the parameters are 
explained by Stephens [4] and by Tan and Ng [23] . 

Since there is no direct experimental data on antiproton energy distribution 
after collision with the air target, we examined several forms of &(E, E') in 
this study. One was very similar to those adopted by Bowen and Moats [1], 
Stephens [4] and Tan and Ng [23]. In this model, Q(E, E') was assumed to be of 
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the form of 1/E'(0 < E < E'), which means the probability that an antiproton 
with initial energy E' possesses energy E after a collision is uniform from E = 
to E — E'. The average energy after a collision is half of the initial energy. 
We call this model a "box- approximation" , because the energy distribution of 
produced tertiary antiprotons has a box-shaped spectrum. Another model was 
proposed by Huang on a different assumption [5]. In his evaluation of total 
inelastic interaction, only annihilation channels were taken into account in the 
inelastic interactions. Non-annihilating inelastic processes were not included 
as a process of antiproton energy loss. This assumption corresponds to taking 
E') = 5(E' — E). This is an extreme case in which a tertiary antiproton 
does not loose its energy in a collision with an atmospheric nucleus. We call 
this model a "^-approximation" . Other intermediate models between the above 
two models could be considered, such as ®(E,E') = l/(E' - E")(E" < E < 
£"), where E" is taken to be various between and E'. The box- and 5- 
approximations correspond to taking E" as and E', respectively. 

In order to compare the observed flux with the calculated results, the follow- 
ing two effects were taken into account: the dependence of the geometrical 
acceptance and antiproton flux on the zenith angle, and the effective observa- 
tion time at each atmospheric depth. The weighted averaged flux taking into 
account these points, (J P (E)), was defined as 



_ J +1 d(cos 9) J* 2 dxJ p (E, x*(x, 9))AT(x u x 2 )ASn(E, cos 9) 



J +1 d(cos 9) \H dxAT(x 1: x 2 )ASn(E, cos 9) 



where J p is a solution of Eq. (1) and x* denotes the effective atmospheric depth, 
which depends on zenith angle of incoming parent particles. AT(x-i,X2) is the 
live observation time between atmospheric depth x\ and x 2 - ASQ(E,cos9) is 
the geometrical acceptance corresponding to the energy and range of zenith 
angle, obtained by the MC simulation. In the calculation for balloon altitude, 
the zenith angle dependence of antiprotons was assumed to be x*(x, 9) = 
xj cos#, because the number of antiprotons produced should be proportional 
to the path length of the primary parent particle. On the other hand, since the 
tertiary antiprotons dominate in those observed at large atmospheric depth, 
they might have collided with nuclei several times and changed their directions. 
In this case, their flux is very sensitive to the mfp inside the atmosphere, 
and their real path length cannot be inferred from the observed zenith angle. 
We simply assumed x* (x, 9) = x in this study. This assumption does not 
significantly change the antiproton spectral shape, but may change its absolute 
flux at a large atmospheric depth. 



6 



6 Results and discussions 



The resultant energy spectra of antiprotons observed at sea level (BESS-1997) 
and at balloon altitude (BESS-2001) are shown in Fig. 3. These results are 
compared with our previous data observed at Mt. Norikura, Japan [7], where 
the mean atmospheric depth was 742 g/cm 2 , together with the results of the 
corresponding model calculations. Since the primary galactic antiproton flux 
was taken into account in the calculation, a sharp edge is seen at the cut- 
off energy of around 3.5 GeV in the calculated flux at 4.5 - 26 g/cm 2 at 
Ft. Sumner. Here, we adopted the galactic antiproton spectrum calculated 
by Mitsui [3] with the solar activity characterized by $ = 1000 MV. The 
antiprotons observed at balloon altitudes are assured to be purely atmospheric 
because the highest energy of the measurement is below the geomagnetic cut- 
off energy. 

The antiproton energy spectrum at small atmospheric depth is mainly deter- 
mined by the production rate of secondary antiprotons, J2a Qa{Ea, x, E) in 
Eq. (1), and the tertiary antiprotons are not dominant. The model calcula- 
tions based both on box- and ^-approximations reproduced the energy spec- 
trum observed at balloon altitude, where secondary antiprotons are dominant 
over tertiary ones. It suggests the production energy spectrum of secondary 
antiprotons was properly treated in our calculation. 

At large atmospheric depth, tertiary antiprotons dominate over secondary 
ones. The calculation results with the box-approximation show good agree- 
ment with observed data at mountain altitude and sea level above 1 GeV. 
The box-approximation was developed by comparing calculated and observed 
proton spectra at sea level and mountain altitude above 1 GeV in previous 
works by Bowen and Moats [1], and Stephens [4]. The observed spectra were 
reproduced well by the box-approximation above 1 GeV, but below 1 GeV the 
^-approximation matches the data more closely. Calculation results with the 
box-approximation show that the energy spectrum is almost flat in a lower 
energy region both at sea level and mountain altitude. However, the energy 
spectra obtained in the ^-approximation decrease below 1 GeV irrespective 
of observation altitude. A model calculation based on a combination of the 
box- and ^-approximations may therefore reproduce the observed spectra bet- 
ter over a wider energy range than the simple box- or ^-approximation. Some 
accelerator experimental results indicate that the shape of the probability 
function <&(E, E') changes depending on the initial energy of a projectile par- 
ticle [3]. Measurement of the atmospheric antiproton spectrum with better 
statistical accuracy over a wider energy range would help to estimate a proper 
shape of $(E,E r ). 
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7 Conclusion 



We measured the atmospheric antiproton spectrum at 4.5 - 26 g/cm 2 and at 
sea level in the kinetic energy range 0.2 - 3.4 GeV for the first time. 

We also referred to our previous observed data at Mt. Norikura in 1999 [7] to 
study the propagation of antiprotons in the atmosphere. The energy spectra 
of antiprotons were calculated for the balloon altitude, mountain altitude and 
sea level. The model calculations based both on box- and ^-approximations re- 
produced the energy spectrum observed at balloon altitude, where secondary 
antiprotons are dominant over tertiary ones. This suggests the production 
energy spectrum of secondary antiprotons was properly treated in our calcu- 
lation. The spectral shapes of our three measurements below 1 GeV were re- 
produced by the ^-approximation, while the calculated flux amplitude at larger 
depths and higher energies was not well-matched to the data. The opposite is 
true for the box-approximation, which matched the flux data above 1 GeV, 
but was not well-matched to the data at larger depths and lower energies. A 
model calculation based on a combination of the box- and ^-approximations 
may reproduce the observed spectra better over a wider energy range than the 
simple box- or ^-approximation. 

Our measurement of antiproton spectra in the atmosphere suggests that the 
shape of the probability function <&(E,E') depends on the initial energy of 
the projectile particle. Measurement of the atmospheric antiproton spectrum 
with better statistical accuracy over a wider energy range is highly desirable 
to improve accuracy of the model calculations. 
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Fig. 1. Atmospheric depth during the BESS-1997 ground observations and the 
BESS-2001 balloon flight. 
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Fig. 2. The identification plots of antiproton events for (a) the BESS-1997 observa- 
tions at sea level, and (b) the BESS-2001 balloon flight. The dotted curves define 
the antiproton mass bands. 
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Fig. 3. The observed antiproton flux at atmospheric depth of 4.5 - 26 g/cm 2 at Ft. 
Sumner in 2001 and at 994 g/cm 2 at Tsukuba in 1997. The flux at 742 g/cm 2 at 
Mt. Norikura [7] is also shown. The results are compared with calculations assuming 
box-approximation (solid lines) and ^-approximation (dashed lines). 
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Fig. 4. Mean free paths assumed in the calculations. 
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Table 1 

Observed antiproton flux at sea level, 994 g/cm 2 at Tsukuba in 1997. 



Kinetic energy (GeV) 
range mean 


Numbers 
antiproton background 


Antiproton flux 
(m~ 2 gr -i s -i QeV -1 ) 


0.20-0.41 
0.41-0.82 0.50 
0.82-1.67 1.19 
1.67-3.40 2.47 


0.0 
2 0.0 
13 0.1 
10 1.9 


8.51 x 10~ 5 upper limit 

6 n2 +6 - 77+0 - 10 x 10~ 5 
D - uz -3.79-0.11 x iu 

n ln +0.70+0.07 1f) -4 
z - iu -0.60-0.07 x iu 

i i i +0.52+0.12 1f) -4 
1-11 -0.44-0.08 x iu 



Table 2 

Observed antiproton flux at the atmospheric depth of 4.5 - 26 g/cm 2 at Ft. Sumner 
in 2001. 



Kinetic energy (GeV) 
range mean 


Numbers 
antiproton background 


Antiproton flux 

(m- 2 sr- 1 s- 1 GeV- 1 ) 


0.20-0.35 0.20 
0.35-0.62 0.54 
0.62-1.09 0.89 
1.09-1.93 1.53 
1.93-3.40 2.61 


1 0.0 
6 0.0 
22 0.0 
50 0.4 
77 4.0 


i c;c; +2.71+0.06 ]r 3 
L - 0D -0.97-0.06 x iu 

4 in +2.29+0.04 -, n -3 
4 - iy -1.52-0.04 x iu 

oa +1.78 +0.24 1f) -3 
O" 34 -1.78 - 0.24 x iu 

1 o fi +0.18+0.05 ln -2 
1 - zu -0.18-0.05 x 1U 

, oi +0.15+0.06 1f) -2 
l -^ L -0.15-0.06 x iu 



14 



